www.gusucode.com > Matlab在化学工程中的应用 > Matlab在化学工程中的应用/实用化工计算机模拟-Matlab在化学工程中的应用/Examples/Chapter 8/OrthExpAnalysis2.m
function OrthExpAnalysis2 % 正交试验的极差分析Variance Analysis of Orthogonal experiment % % Author: HUANG Huajiang % Copyright 2003 UNILAB Research Center, % East China University of Science and Technology, Shanghai, PRC % $Revision: 1.0 $ $Date: 2003/07/12 $ % % [Ref] 何为和陈际达编,优化试验设计法及其在化学中的应用,成都: % 电子科技大学出版社,1994(p.20) NL = 2; % 水平数 NF = 6; % 因素数 r = 4; % 每个水平有r次试验 DATA = ... % DATA矩阵的前三列是正交试验方案,第四列为试验指标--收率 [1 1 1 1 1 1 65; 1 1 1 2 2 2 74; 1 2 2 1 1 2 71; 1 2 2 2 2 1 73; 2 1 2 1 2 2 70; 2 1 2 2 1 1 73; 2 2 1 1 2 1 62; 2 2 1 2 1 2 67]; A = DATA(:,1:NF); Y = DATA(:,NF+1); for i=1:NL % i:水平 for j=1:NF % j:因素(列号) K(i,j) = sum(Y(find(A(:,j)==i))) end end k = K/r % 极差 R = max(k(:,:)) - min(k(:,:))